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I. INTRODUCTION 



The Space Shuttle, operated by the National Aeronautics and Space Administration 
(NASA), is the only manned space flight vehicle currently operated by the United States. The 
Space Shuttle consists of the Orbiter, External Tank (ET), and two Solid Rocket Boosters 
(SRB’s). The Orbiter contains the astronauts and payloads. It is the only piece of the Space 
Shuttle that completes the entire mission. The Orbiter returns from space to land similar to a 
conventional gliding airplane. The continued safe operation of the Space Shuttle is the number 
one concern of the Space Shuttle Program (SSP). A recent safety concern is the directional 
control of the Orbiter during one segment of the landing/rollout phase [Ref. 1]. This perceived 
directional control deficiency is described below. 

A. BACKGROUND 

During the landing and rollout phase of flight, the gains associated with the piloting task 
are greatly increased. The margin for error is vastly decreased as pitch, roll, and yaw requirements 
become more stringent through the flare-to-touchdown and subsequent rollout. During this 
critical phase, pilot control commands will be more frequent and often of greater magnitude to 
make timely corrections to flight path and counter external perturbations such as crosswinds, wind 
shear, or gusty conditions. These increased demands on the control system can occasionally result 
in a condition known as Pilot Induced Oscillation (PIO). 

PIO occurs when a pilot attempts a particular control task which is beyond the capabilities 
of the aircraft control system. PIO can only occur with the pilot actively in the vehicle control 
loop and is therefore not a stability concern. A PIO is typically manifested by an uncommanded 
oscillatory motion that can only be stopped when the pilot reduces control inputs or ‘backs out of 
the loop.’ The severity of PIO can range from an annoyance to complete loss of control. The 
result of a PIO situation can often be unpredictable. What might be an annoyance during landing 
one day, might be disastrous on the next. PIO can result from many different root causes 
including an excess amount of friction in a mechanical control system, actuator rate limiting, or 
transport delays in today’s most advanced flight control systems. In all cases, a large phase lag 
from the control input to aircraft response, lies at the heart of the problem. 

The Orbiter incorporates a ‘fly-by-wire’ control system to command the flight control 
surfaces, depicted in Figure 1. 
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Inboard Eievon 




The flight controls, control stick, rudder pedals, and speed brake switch are not 
mechanically connected to their respective control surfaces. Orbiter control is afforded through a 
myriad of feedback control loops, filters, gain multipliers, integrators, differentiators, etc. The 
goal of these components, in conjunction with a feel system, is to supply the pilot with an Orbiter 
which responds acceptably to all control inputs. Fly-by-wire control systems, if not carefully 
designed, can be susceptible to PIO. 

The directional control problem of interest is most evident after touchdown while the 
Orbiter is in the two point stance with the mainwheels on the ground and the nosewheel in the air. 
This flight segment, from initial touchdown until derotation is complete, is approximately 20 
seconds in duration. With dispersed touchdown speeds between 185 and 225 knots, a small 
deviation in track can result in a runway departure very quickly. Recent Orbiter landings have 
shown that there is a directional PIO problem in the two point stance[Ref. 1], The problem seems 
to abate after nose gear touchdown due to the additional directional control afforded by the Nose 
Wheel Steering (NWS) system. The PIO is associated with the task of capturing and tracking the 
runway centerline and may be exaggerated in conditions of greater crosswind. The possibility of 
the Orbiter becoming directionally uncontrollable and departing the runway is unacceptable. This 
problem raises some critical questions including: 
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1 . Are there operational impacts? Must the crosswind limit be reduced to protect against 
PIO? Should Shuttle landings require wider runways and/or lakebed landings? 

2. Are there hardware modifications such as changes to the rudder actuators that would 
improve directional control? 

3. Are there potential software modifications to the control laws which will reduce or 
eliminate this PIO condition? 

4. If PIO is evident in a nominal landing, what would happen in the case of a blown 
main gear tire? 

1. Operational Impacts 

Reducing the crosswind limit may reduce the potential for direction PIO; however, it will 
not eliminate the problem. In fact, the landing which experienced the most significant directional 
PIO had a 2 knot crosswind [Ref. 2], Other flight data also shows that even with minimal 
crosswind, directional PIO can occur. 

In addition, there is appreciable pressure to increase the crosswind limits for the following 
reason. The commitment of the United States to building the International Space Station has 
extended the mission of the SSP from pure science to upholding national political obligation. The 
proposed International Space Station will be in a high inclination orbit which will reduce launch 
opportunities and decrease the length of launch windows. Landing weather conditions are a 
concern on launch day since the potential for an emergency landing such as a Return To Launch 
Site (RTLS) or Trans-Atlantic Landing (TAL) can occur on launch day. Meeting the Space 
Station construction schedule would be much easier if launches and landings could be performed 
in a greater variety of wind conditions. 

Requiring lakebed landings for all Shuttle flights is unacceptable. Orbiter landing mass 
properties are greatly restricted for lakebed landings to maintain nose gear loads within limits 
[Ref. 3], Most Orbiter landing mass properties require concrete runways. Further, landing on a 
lakebed requires landing at either Edwards Air Force Base (EAFB) in California or Northrup, 
New Mexico which is impossible for most ascent aborts modes. 

The Shuttle Landing Facility (SLF) at Kennedy Space Center (KSC) and the concrete 
runway used for Shuttle landings at EAFB are 300 ft. wide. The SLF has 50 ft. of usable shoulder 
on each side. In comparison, some contingency and abort runways are less than 150 ft. wide. 
Obviously, it is impossible for the SSP to widen all runways available for landing. Restricting the 
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Orbiter from landing on runways considered too narrow, greatly reduces abort coverage (no where 
to landing with an SSME failed) to an unacceptable level. 

2. Hardware Modifications 

Hardware modifications, such as redesign of the rudder actuator to increase deflection 
rate, may reduce susceptibility to PIO. However, hardware modifications are expensive, require 
recertification, and demand vehicle down time to implement the changes. 

3. Software Modifications 

The obvious long term solution to the PIO problem lies in redesigning the existing flight 
control laws to afford better directional control in a wider variety of wind conditions. 

A redesign of this magnitude would traditionally require the control engineer to estimate 
the new control laws followed by experimental Verification and Validation (V&V) of the new 
control laws in the Vertical Motion Simulator (VMS) located at the NASA Ames Research Center 
(ARC) located in Sunnyvale, California. This procedure typically entails qualified astronauts 
traveling from the Johnson Space Center (JSC) in Houston, Texas to the ARC VMS to fly 
numerous simulated approaches in varying wind conditions. This evolution can be time 
consuming and costly. 

With the advent of faster more powerful personal computers and workstations, and the 
availability of very capable software packages which can operate in both the frequency and time 
domains, the possibility of writing and evaluating new control laws immediately has become a 
reality. This breakthrough can eliminate the costly iterative process of tweaking and reevaluating 
control laws in an expensive motion base simulator. The final V&V is still necessary, but the total 
time and cost are drastically reduced. 

The use of simulation to evaluate control laws rely heavily on the numerous models used 
to develop the simulation. A normal flight test program would not rely solely on simulation 
results for solutions to PIO problems. The Shuttle Program must rely on these simulations results 
due to the nature, cost, and risk of flying Shuttle flights only to evaluate control law changes. 
Flight data is gathered from each flight; however, flight conditions do not test the extreme ends of 
the envelope. There is no possibility of performing a complete flight test program. 



4 



B. 



PURPOSE 



This thesis has two objectives: 

(1) Develop changes to the control laws to reduce the directional PIO tendency present 
between mainwheel touchdown and completion of derotation while the Orbiter is in 
the two point attitude. 

(2) Evaluate these changes using a Commercial Off-The-Shelf (COTS) PC-based 
software package. 
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II. ORBITER DIRECTIONAL STABILITY AND CONTROL 



The nominal Orbiter rollout profile, directional stability, and directional control in the two 
point stance are described in this chapter. 

A. NOMINAL ROLLOUT PROFILE 

The nominal touchdown and rollout profile, with approximate speeds and time from 
touchdown, is described in Table 1 [Ref. 2, 3 and 4], 

Table 1. Nominal Touchdown/Rollout Profile 



Event 


Airspeed (KEAS) 


Time (sec.) 


Touchdown 0 * 


205 (195) 


0 


Brakes functional^ 1 


200 


1.9 


Drag chute deploy* 3 ' 


195 


3.0 


Commence auto derotation <4) 


185 


11.0 


Commence manual 
derotarion (5> 


175 


16.0 


Drag chute fully disreefed 


175 


16.0 


Nosewheel touchdown 


145 to 155 


19.0 


First brake application* 6 * 


140 (KGS) 


23.0 


Drag chute Jettison* 7 ' 


60 ±20 (KGS) 


34.0 


Wheels stop 


0 


50.0 









Notes: 1 . Targeted touchdown speed of 205 KEAS for Orbiter gross weights > 

220K lbs. (195 KEAS for Orbiter weights < 220K lbs.). 

2. Built In Test (BIT) requires 1.9 sec. after Weight On Wheels Lock-On 
(WOWLON). 

3. Chute deploy after main gear touchdown. 

4. Nominal derotation is automatic with beep trim activation targering 1.5 
deg/sec. 

5. Only performed if beep trim fails. 

6. Brake initiation at 140 KGS or 5,000 ft. runway remaining, whichever occurs 
first. 
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7. If drag chute is not jettisoned by 40 KGS, the chute will remain on to 
minimize damage to main engine bells. 

B. DIRECTIONAL STABILITY IN THE TWO POINT ATTITUDE 

Directional stability in the two point stance is different than the directional stability of the 
Orbiter in flight. In flight, the Orbiter would typically rotate about the Center Of Gravity (COG). 
In the two point stance, assuming no lateral skidding of the mainwheel tires, the Orbiter is 
constrained to rotate directionally about a point on a line that runs laterally through the areas of 
contact between the mainwheels and the runway. This point is referred to as the Center Of 
Rotation (COR). Directional stability is governed by four main factors in the two point attitude: 

(1) The relative position of the COG to the COR. 

(2) The increased friction on the downwind mainwheel caused by the rolling moment 
generated as a result of the dihedral effect or roll due to sideslip, Lp, generated 
primarily by the highly swept wing. 

(3) The aerodynamic ‘weathercock’ effect or yaw due to sideslip, Np, generated primarily 
by the large vertical tail. 

(4) The relative position of the drag chute to the COR. 

The first factor will attempt to keep the Orbiter traveling in a straight line with respect to 
the runway in a zero crosswind condition but can act to direct the Orbiter away from a crosswind 
as a result of the second factor. The third factor will tend to fair the Orbiter into a crosswind in 
flight. After touchdown, with the COR shifting aft, this ‘weathercock’ effect becomes a much 
weaker contributor. The forth factor will attempt to fair the Orbiter into a crosswind. 

1. Center of Gravity (COG) Effects 

The nominal COG envelope is 32 in. long and 3 in. wide, centered laterally and 
positioned approximately 8 !4 ft. forward of the mainwheels [Ref. 5]. The COG will always be 
well ahead of the mainwheels in the two point stance. The COG geometry with respect to the 
mainwheels is depicted in Figure 2. 
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With this geometry, the COG acts as a pendulum in the horizontal plane. The 
deceleration of the Orbiter caused by the friction in the mainwheels and the mainwheel brakes acts 
on the COG to produce a force in the Orbiter’s direction of motion. This force is analogous to the 
gravitational force on a normal, vertically oriented, pendulum. The result is a restoring yaw 
moment which is always attempting to maintain the Orbiter on a straight track in the direction 
from the COR through the COG. In a zero crosswind condition with no mainwheel braking, the 
COR will most likely be positioned directly between the mainwheels. In this event, the Orbiter 
will have positive directional stability in the direction of motion. In a crosswind case, the high 
dihedral of the Orbiter can become a directionally destabilizing factor. 

2. Dihedral Effect 

The Orbiter has a large dihedral effect due primarily to the highly swept wing. Dihedral is 
the tendency for an aircraft to roll away from sideslip. In a landing scenario, the sideslip angle can 
be calculated as follows: 
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crosswind component 



Orb iter Indicated A ir speed 



The rolling moment generated in a highly swept wing aircraft in response to sideslip can 
be explained by two factors. First, the upwind wing senses flow closer to perpendicular to the 
leading edge which decreases flow in the spanwise direction; thereby increasing lift. The 
downwind wing has the opposite effect of increased spanwise flow; therefore decreased lift. The 
second, compounding effect, is a blanking effect on the downwind wing caused by the large 
fuselage structure of the Orbiter. This reduction in airflow also decreases lift. 

Another contributor to positive dihedral is the high vertical tail. The center of pressure 
associated with the airflow across the vertical tail is well above the COR which would cause a 
further roll away from crosswind. The net result is a rolling moment away from crosswind as 
depicted in Figure 3. 



This rolling moment will cause an increased down force on the downwind strut and a 
lighter force on the upwind strut. This differential load will cause increased friction on the 
downwind wheel and a subsequent shift in the directional COR to the downwind side of Orbiter 




High Vertical Tail 



Side Force from 



Figure 3. Orbiter Dihedral 
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centerline. This shift in the COR, in conjunction with the COG effect discussed in subsection 1, 
will act to cause the Orbiter to deviate away from crosswind as depicted in Figure 4. 




3. Weathercock Effect 

The weathercock effect is typically a tradeoff between the destabilizing effect of the 
fuselage and the stabilizing effect of the vertical tail. The aft location of the COG of the Orbiter 
in flight is somewhat destabilizing since the stabilizing moment arm associated with the fin is 
shorter and the destabilizing moment arm of the fuselage is longer. When the Orbiter attains the 
two point stance, the COR shifts even further aft. This change creates a condition of weak 
instability. The Orbiter will tend to point away from a crosswind as depicted in Figure 5. 
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Figure 5. Orbiter Yaw Due to Weathercock 



4. Drag Chute Effect 

The drag chute was incorporated on to the Orbiter to decrease landing rollout distances for 
certain contingency runways. Along with the decreased rollout, there is a small amount of 
positive weathercock effect. It is important to note that the drag chute is nominally fully disreefed 
three seconds prior to nose wheel touchdown as described in section II. A [Ref 2 and 3]. 
Therefore, the drag chute will only be a brief factor to steady-state directional stability in the two 
point stance. Yet the three second impulse may be enough to perturb the steady case. The drag 
chute will blow downwind from a crosswind. Since the drag chute is attached aft of the center of 
rotation, the resulting yaw moment on the Orbiter is into the crosswind as depicted in Figure 6. 
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5. Frequency and Damping 

A summary of the directional stability parameters is shown in Table 2 along with their 
relative strengths. 



Table 2. Summary of Directional Stability Parameters 



Parameter 


Type 


Relative Strength 


Direction 


COG 


Inertial 


Strong 


In direction of motion 


Dihedral 


Aerodynamic 


Weak 


Away from crosswind 


Weathercock 


Aerodynamic 


Weak 


Away from crosswind 


Drag Chute 


Aerodynamic 


Weak 


Into crosswind 



All of these factors combine to yield fairly weak directional stability in the two point 
stance. This weak stability, when modeled with a simple second order system, would yield a low 
open loop frequency response to a disturbance. While this data point has never been collected in 
flight, a rough estimate of what one might call ’the directional natural frequency', o>l>> can be made 
by observing video of actual closed loop cases of PIO. The videoThe open loop natural 
frequency, co D , is estimated to be on the order of 1 rad/sec. The only yaw damping, N r , is most 
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likely due to the large vertical tail and possibly any differential friction in the mainwheels. Both 
of these factors would be small which yields very low directional damping. The complete 
dynamic stability of the Orbiter in the two point stance can be modeled by the following second 
order equation where Njd is total damping and Nts is total stability. 

s 2 + N m s + Nj S 

where co D = yj N IS 

Njd 

and the damping ratio, <% D = — 

2g> d 

The natural frequency and damping ratio will be estimated from the simulation results. A 
low frequency and low damping ratio will indicate the tendency of the FCS to be susceptible to 
PIO. 

C. DIRECTIONAL CONTROL IN THE TWO POINT ATTITUDE 

Directional control in the two point stance can be affected with two controllers; 
aerodynamically with the rudder, and mechanically with differential braking. Brakes are not 
employed until derotation is complete to avoid excessive nosewheel slapdown rates and to reduce 
the risk of a blown main tire. Braking nominally begins at 140 knots [Ref. 3], Differential 
braking would only be employed in extreme situations, such as a flat main gear tire. Therefore, 
directional control in the two point stance would be afforded only by the rudder. 

1. Rudder 

The rudder is actually a dual purpose control surface acting as both a directional controller 
and a speedbrake as depicted in Figure 7. 
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Figure 7. Orbiter Rudder/Speedbrake 



The surface is mechanically limited to a maximum deflection rate of 16 deg/sec normally 
[Ref. 6]. Postlanding the deflection rate is Anther limited tol3deg/sec since the speedbrake is full 
open. The software limits the deflection rate to 14 deg/sec so that sufficient hydraulic fluid is 
available to other control surfaces. In the event of the lost of one or two hydraulic system, the rate 
is limited to 12 deg/sec and 7 deg/sec, respectively [Ref 4]. The surface areas is approximately 
97 ft. 2 and can generate up to 2.25 lateral g's. 

2. Directional FCS Susceptibility to PIO 

The FCS has many limitations which result in the high potential for directional PIO in the 
two point stance. The slow surface rate of the rudder, slow response time, high control power of 
the rudder due to the large surface area, and rudder placement with respect to the gear are 
conducive to PIO. If a large pedal input is sustained until the yaw rate is well underway, loss of 
control is probable. The slow surface rate does not allow an opposite input to be effected quicky 
enough to counteract the yaw rate. 

In addition, ARC VMS studies have indicated rudder rate saturation is encountered in 
high crosswind landings or in performing postlanding lateral maneuvers stress tasks [Ref. 7]. 
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III. RUDDER CONTROL 



This chapter describes the rudder control path including the software command path and 
rudder actuator. Detailed drawings and description of the Orbiter FCS control laws are found in 
the Software Requirements Document [Ref 8] and Lockheed drawings [Ref. 9]. 

A. RUDDER COMMAND PATH 

The control logic for the rudder surface may be broken down into two basic steps. The 
first converts rudder pedal deflection angle into lateral acceleration error command. The Orbiter’ s 
flight control system uses lateral acceleration feedback to null errors to the commanded lateral 
acceleration and provide directional stability augmentation. The second step converts the lateral 
acceleration error command into a yaw rate which is filtered of Orbiter structural bending modes 
and finally converted into a command to the control surface for degrees of deflection. This 
command is sent to Priority Rate Limiting (PRL), which limits the rudder deflection rate 
depending on the number of hydraulic systems available. These two steps will be henceforth 
referred to as the rudder command feedback and rudder deflection command portions. 

1. Rudder Command Feedback 



A simplified block diagram of rudder pedal position conversion to lateral acceleration 
error command is shown in Figure 8. 
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Figure 8. Rudder Command Feedback 
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The various components of the rudder command feedback path are discussed in the 
subsections below. 

a. Rudder Pedals 

There is a Rudder Pedal Transducer Assembly (RPTA) located at both the 
Commander and Pilot stations. Each RPTA consists of two pedals and a transducer which 
converts degrees of rudder pedal deflection into electronic signal. The RPTA uses a time constant 
of 100 ms in the mechanical to electronic transfer. This signal is processed at 6.25 Hz and passed 
to the rudder command shaper with the name tag DRMAN. 

b. Rudder Command Shaper 

The rudder command shaper converts DRMAN into a shaped and scaled rudder 



pedal command named DRMS as shown in Figure 9. 




Figure 9. Rudder Command Shaper 



The shaping schedule is parabolic and conforms to the following function: 
DRMS = ( 0.042)DRMAN 2 + (Q.l2l)DRA4AN . A maximum deflection limit of 22.5 degrees 
is imposed on the output signal. 

c. Rudder Pedal Gain Scheduler 

The rudder pedal gain scheduler produces a gain signal, which when multiplied 
with the shaped and scaled rudder pedal command, DRMS, produces a lateral acceleration 
command. This gain is scheduled as a function of True Airspeed (TAS) and is given the name 
GNYDRM as depicted in Figure 10. This gain also converts the units of the signal from degrees 
of pedal deflection to lateral g’s. The rudder pedal gain scheduler samples at 1 .04 Hz. 
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Figure 1 0. Rudder Pedal Gain Scheduler 



d. Rudder Trim Command 

The rudder trim command signal is generated by converting the command signal 
from the rudder trim panel into a commanded lateral acceleration bias. This is performed with a 
rectangular integrator which ramps the output signal and also limits the maximum bias to 0. 1 g's. 
This signal is sampled at 6.25 Hz. The output signal is the lateral acceleration trim signal, DRTI, 
which is summed with GNYDRM from the rudder pedals to produce a total lateral acceleration 
command, DRTMS, as depicted in Figure 11. 



GNYDRM x DRMS 
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e. 



Lateral Accelerometers 



There are four Accelerometer Assemblies (AA) located forward of the crew 
compartment. These AA’s directly measure lateral acceleration and simultaneously filter the 
signal as shown below: 

Measured Acceleration 



15 2 



\t 



146 2 



— » Lateral Acceleration Signal (NYP) 



s 2 +2(.5)(15)s+15 2 As 2 +2(.8)(146)s+146 2 



This signal is sampled at 25 Hz and passed to the lateral acceleration feedback 

compensator. 

f Lateral Acceleration Feedback Compensator 

The lateral acceleration feedback compensator compares the actual lateral 
acceleration, NYP, as passed from the AA’s, to the commanded lateral acceleration, DRTMS, to 
generate an error signal, or additional lateral acceleration command. A low pass filter is 
incorporated in this procedure to smooth the output. The resulting signal is sampled at 25 Hz. 
The process is depicted in Figure 12. 
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Figure 12. Lateral Acceleration Feedback Compensator 



2. Rudder Deflection Command 

A simplified block diagram of lateral acceleration conversion to rudder deflection 
command is shown in Figure 13. 
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Figure 13. Rudder Deflection Command 



The various components of the rudder deflection command path are discussed in the 
subsections below. 

a. Lateral Acceleration Gain 

The lateral acceleration gain scheduler produces a gain signal, GRAY, which 
when multiplied with the lateral acceleration error command, DRTMS, produces a yaw rate 
command, DAY. This gain is scheduled as a function of Mach as depicted in Figure 14. This 
gain also converts the units of the signal from g’s to deg/sec. The sample rate for this scheduler is 




Figure 14. Lateral Acceleration Gain Schedule 



b. Cosine of the Angle of Attack (a) 

A dimensionless signal, cos(a), is mixed with the yaw rate command signal and 
summed with the yaw rate error to produce a yaw rate command about an inertial axis. This 
signal will always reduce the yaw rate command signal in producing a yaw rate with respect to the 



21 




runway. The cos(a) is sampled at 1 .04 Hz. The output signal from this mix is referred to as the 
unfiltered yaw rate command, DRRCLJF. 

c. Rate Gyro Assemblies 

There are four Rate Gyro Assemblies (RGA’s) located on the aft bulkhead below 
the floor of the payload bay. Each of the RGA’s sense rotation about the three vehicle axes and 
simultaneously filter the signal as shown below: 

Measured Acceleration 

— 2 y - > Body Rates ( P,R,Q ) 

x 2 + 2(.7)(420 ).y + 420 2 J 



50 A 
vx +50 J 



These rates are the primary feedback to the flight control system to maintain 

control. 

d. Yaw Rate Structural Filter 

A second-order structural bending filter is applied to the yaw rate feedback. 

e. Yaw Rate Error Determination 



The filtered body yaw rate, DRPRM, is summed with the product of the sideslip 
angle times the cosine of the angle of attack, BCOSALF, and yaw rate command, DAXFDC, to 
determine the yaw rate error, DRERR, as shown in Figure 15. The sideslip angle is not used in 
the approach and landing flight phase. Therefore, BCOSALF is zero in this analysis. The 
resulting signal is sampled at 25 Hz. 



BCOSALF DAXFDC 





DRPRM 






DRERR 



Figure 15. Yaw Rate Error Determination 
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f Structural Bending Filter 

A structural bending filter is included in series to remove any known Orbiter body 
bending modes. The resulting signal is the yaw command, DRRC, which is processed at 25 llz. 
The filter is made of two consecutive second order notch filters. This fourth order system will 
inherently result in transport delays which would only contribute to the time lag problem. 

g. Yaw Rate Gain 

The yaw rate gain scheduler produces a gain signal, GDRC. This gain and a 
constant, GERYFWDS, is multiplied with the yaw command signal, DRRC, to produce an initial 
rudder deflection command, DRCPF. This gain is scheduled as a function of Mach and dynamic 
pressure as depicted in Figure 1 6. This gain also converts the units of the signal from deg/sec to 
deg. The sample rate for this scheduler is 1 .04 Hz. 




Figure 16. Yaw Rate Gain Scheduler 
fu Rudder Deflection Limiter 

The rudder deflection limiter provides a maximum rudder deflection limit in the 
positive and negative directions. The deflection limit, DRC LIM, is 24. 1° and 27. 1° in flight and 
post Weight On Wheels Lock-On (WOWLON), respectively. Within these limits, the limiter 
passes a one-to-one correlation from the filtered rudder deflection command, DRCPF, to the 
actual rudder deflection command, DRC. The WOWLON signal is sampled at 6.25 Hz, but the 
rudder deflection limiter processes at 25 Hz. The system is configured as shown in Figure 1 7. 
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DRCPF 




Figure 1 7. Rudder Deflection Limiter 
/. Priority Rate Limiting (PRL) 

The rudder deflection command is sent to PRL. PRL limits control surface 
deflection rates to manage the amount of hydraulic fluid sent to each of the actuators. There are 
three Auxiliary Power Unit (APU)/hydraulic systems available to provide hydraulic fluid to the 
control surface actuators. Normally, with three hydraulic systems available, PRL limits the rudder 
deflection rate to 14 deg/sec [Ref. 4], This slow deflection rate is conducive to the pilot’s rudder 
inputs being out of phase with Orbiter response. The deflection rate is further reduced depending 
upon the number of hydraulic systems failed. With one and two systems failed, the rudder 
deflection rate is limited to 12 and 7 deg/sec, respectively. 

B. RUDDER ACTUATOR 

The rudder/speedbrake actuators receive hydraulic fluid from all three hydraulic systems 
through a switching valve. The switching valve produces one hydraulic output pressure source for 
the rudder/speedbrake servovalves. The power drive unit, mechanical actuator, and actuator drive 
shaft are shown in Figure 1 8. 
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Figure 18. Rudder/Speedbrake Actuator Location [Ref. 6] 
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VERTICAL STABILIZER 
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IV. EQUATIONS OF MOTION AND STABILITY AND CONTROL 



DERIVATIVES 



A. ASSUMPTIONS 

The following assumptions were used in developing the Equations of Motion (EOM) for 
the Orbiter in the two point stance: 

1 . Inertia] frame fixed to the flat Earth since the Orbiter is on the ground. 

2. Constant vehicle mass properties since there is no propellant consumption and 
nothing ejected until drag chute deploy. 

3. Vehicle is symmetric about the XZ plane. Therefore, 1^ = 1^ = Cyo = Cio = C no = 0. 

4. Equations of motion are linear about a steady state flight condition. Small 
perturbation theory implies that the squares and products of perturbations may be 
neglected. For every perturbation angle a. Sin (a) = a and Cos (a) = 1 . 

5. Atmospheric property variations are negligible. Air density, temperature, pressure 
changes are negligible over a twenty second period and less than 10,000 feet rolling 
distance on the ground. 

6. Steady state rotation rates, side velocity, and bank angle are all zero. The Orbiter in 
the two point stance will not perform rapid maneuvers or fly with significant steady 
state sideslip. 

7. The air flow is quasi-steady. In other words, the pressure distributions adjust instantly 
to rate changes. This assumption is generally valid for angle of attack rate a if the 
following condition is satisfied [Ref. 10]: a c/2U| < 0.04. For the Orbiter, the chord 
length c = 39.56 ft. and the initial velocity Ui = 353 ft./sec. The inequality is 
satisfied if a < 0.714 rad/s = 40.9 deg/s. In the two point stance prior to the initiation 
of derotation, a is essential constant. After initiation of derotation, a is 1.5 deg/sec. 
The stability derivatives with respect to a may be neglected. 

8. Motion in the Z direction is constrained due to the runway. Therefore, the downward 
velocity W=0 and the associate forces and moments are zero. 

9. Vehicle is at constant a. Therefore, all pitching moment derivatives are zero, the 
pitch rate, q, is zero, and all control derivatives with respect to q are zero. 
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10. The elevons are at trail (0°) at touchdown and remain there until the initiation of 
derotation [Ref. 2], Therefore, all variations of the coefficients due to elevator and 
aileron are zero. 

1 1 . Forward velocity, U is constant. 

B. EQUATIONS OF MOTION 

The general EOM in the body-fixed coordinate system obtained from Roskam [Ref. 10] 
are defined below. The subscript denotes the axis and the apostrophe denotes time derivative. All 
symbols are defined in the List of Symbols, Acronyms, and/or Abbreviations on page xii. 
Although the control problem under investigation is a directional control problem, both the 
longitudinal and lateral/directional equations are included in the derivation and simulation for 
completeness. The general force, moment, and kinematic equations are: 

F x = m(U -VR+WQ)+mgSin(0) 

F y = m(V +UR-WP)-mgCos(©)Sin(0) 

F z = m(W -UQ+VP)-mgCos(<l>)Cos(0) 

L = I X xP'-IxzR'-IxzPQ+(Izz-Iyy)RQ 
M = I yy Q+(I>orWPR+UP 2 -R 2 ) 

N = IzzR-IxzP +(Iyy-Ixx)PQ+IxzQR 
P = 0 -vp Sin(0) 

Q = 0 Cos(0)+vp Cos(0)Sin(0) 

R = vp Cos(0)Cos(O)-0 Sin(O) 

Each variable is the sum of a steady state components, denoted with a subscript “1”, and 
perturbation components, denoted with lower case: 



U = Uj+u 


V = V,+v 


W= Wj+w 


P = Pi+P 


Q = Q]+q 


R = R,+r 


4 / = 4 > ,+vp 


0 = ©,+e 


0 = 0|+<j) 


v P’ = v P,’+vp- 


©' = ©,'+0’ 


0 = 0] +<j) 


Fx + F xJ + f x 


Fy + FyJ + fy 


F z + F z) + f z 


L = L,+l 


M = M ) rtn a 


N=Nj+n 



The steady state forces, moment, and kinematic equations are as follows: 
F x i = m(U i-ViR|+WiQj)+mgSin(©i) 
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F y i = m(V i+UiR|-W 1 Pi)-mgCos(0i)Sin(d>i) 

F 2| = m(W|-U iQi+V ,Pi )-mgCos(<t>, )Cos(0 1 ) 

Ll = IxxP 1-IxzR rl-v-P I Ql +(Izz-l>y)R|Ql 

M, = lvvQ', +(Ixx-lzz)PlR|+lxz(Pl 2 -R. 2 ) 

N, = I^R 1 -l^P 1 +(I >y -Ixx)P 1Q1 +IxzQlRl 
P] = d> i -vg ]Sin(0) 

Qi = 0 , iCos((P])+v(/ 1 Cos(0i)Sin(Oi) 

Ri = v|/ !Cos(©i)Cos((|)|)-0 !Sin(Oi) 

The pitch and bank angles are defined such that: 

Sin(0) = 0 Sin(<t>) = <)> Cos(0) = Cos(<(>) £ 1.0 

This small angle assumption restricts the pitch and bank angle perturbations to less than 
15 degrees. This limit is reasonable for landing and rollout and does not constitute any serious 
restriction. Perturbations are assumed to be sufficiently small such that products and cross- 
products of the perturbations are negligible with respect to the perturbation itself: 
pv = pw = pq = pr = qu = qw = qr = ru = rv = 0 
u 2 = v 2 = w 2 = p 2 = r 2 = <p0 = vg0 = v(/(p = (p0 = 0 

The perturbation force, moment, and kinematic equations are as follows: 
f x = m(g0Cos(0)+ W i q- V i r+u -Ri v+Qi w) 

f y = m(-g(pCos(Oi)Cos(0i)+g0Sin(cD 1 )Sin(0i)-W|p+Uir+RiU+v -Piw) 
f z = m(g(pCos(0 1 )Sin(cp|)+g0Cos(d)i)Sm(0|)+V 1 p-Uiq-QiU+PiV+w ) 

1 = Up -(QiP+Piq)Ixz + (-Iyy + lz2)(Riq+Qir)-Ix2r 
Itla — 21xzP 1 P+(Ixx-I«)R|P+lyyq "KIxx"lzz)P lf-H^RiT 

n = (-ixx+4y)(QiP + Piq)-lxzP ^(Riq+QiOixz+lzzt 
p = -04^ i Cos(0 ] ) - vp Sin(0i)+cp 

q = 0 Cos(d>])+vj/Cos(0 1 )Sin(d> 1 )-0 v P ,Sin(d) 1 )Sin(©])+4 > ,Cos(d>,)Cos(0i)-0 iSin(d>,))(p 
r = -(0 |Cos(d>i)+ v P 1 Cos(©i)Sin(Oi))(p+v}/Cos(Oi)Cos(0i)-0 Sin^iJ-OT iCos(<t>i)Sin(0i) 

1. Longitudinal Equations of Motion 

The longitudinal forces in the stability axis system are defined as: 

F X 1 D F zl L Dj — CoiqdynS Lj — CuqdynS — C m iqdynSc 

The non-dimensional stability and control derivatives in the stability axis are defined as: 
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The longitudinal dimensional stability axis derivatives are defined as: 
Axial force: 

-^,^. +2C 0 ft 



x.. = 



( ldyn^{^D u + D , ) 

mil x 



sec 



X„ = 



m 



v =4 1 



rad sec 2 M U x sec 



x ^4^ SCC D, ft _ ^dyn SC D it ft 

q 2mU x rad sec s * 



m rad sec" 



Normal force: 

~ ( i cfyn^{^L, + ^Q, ) 1 ^ - Q dynSjf-' L a + ) ft_ 



Z.. = 



mU x sec 



Z. = 



z„ 1 



m 



rad sec 2 K U x sec 



z _ -<i*, ScC i, fi 

2 mU x rad sec 






ft 

m rad sec" 



Pitching moment: 

+2( -Mi) rad Iy q^ScC^ ] 



M.. = 






ft sec 1 



AT... = 



>y 



sec 



M a rad 
U x ft sec 



Rdytfic C ] (JdyJcC 1 

M, = _ * — M s = 



q 21 iy U x sec 



lyy sec : 



Recall that motion in the Z direction is constrained due to the runway. Therefore, the 
downward velocity W=0 and the associate forces and moments X u , Z u , and M vv are zero and not 
included in the state space equation. 

The pitch moment of inertia changes slightly when the Orbiter touches down. The center 
of rotation shifts from the center of mass to the axis formed by the line between the points where 
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the main gear tires meet the runway. The parallel axis theorem [Ref. 1 1] allows a simple shift of 
the moment of inertia with the following equation: I ymg - 1 yy + mil 2 where l yrng is the moment of 

inertia about the main gear, m is the mass of the Orbiter, and h is the distance between the center 
of mass and the new rotation point. All pitch moment dimensional stability derivative equations 
will be adjusted accordingly by multiplying by I >y /I ymg . 

The following are the perturbed longitudinal equations of motion: 

u'= X u u+X a a + X q q -gcos(0,)#+ X s 8 e 

, Z u u Z a (U\ + z q )q gsinf©,)# Z s 8 e 
(/, U x (/, U y U x 

q'= M u u+ M a a + M q q + M s 8 e 

0=q 



In state space form, X’=AX+BU: 
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2. Lateral-Directional Equations of Motion 



The lateral-directional forces are defined as: 

Fyi Lj Cuq d>1) Sb Ni C n^dynSb 

The non-dimensional stability and control derivatives are defined as: 
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The lateral-directional dimensional stability and control derivatives are defined as: 
Side force: 
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The yaw moment of inertia changes slightly when the Orbiter touches down. The center 
of rotation shifts from the center of mass to a point between the main gear. It is assumed that this 
point will be directly between the main gear. The parallel axis theorem [Ref. 1 1] allows a simple 

shift of the moment of inertia with the following equation: / = /„ + mh 2 where I 7rog is the 

moment of inertia about the main gear, m is the mass of the Orbiter, and h is the distance between 
the center of mass and the new rotation point. All yaw moment dimensional stability derivative 
equations will be adjusted accordingly by multiplying by 

The roll moment of inertia is assumed to remain referenced to the center of mass as the 
main gear struts are free to compress and extend which permits the Orbiter to roll about a point 
near the center of mass. I xx will not be adjusted. 
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The following are the perturbed lateral-directional equations of motion: 



P=Y f 



A 

u, 



+ F 



U, 



+ Y 



1 



r U, 
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For small 0’s and small 1^, inertial coupling is negligible; therefore' the Ixzt’/Ixx and I 0 p’/I vx terms 
are disregarded. 



In state space form, X’=AX+BU: 
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3. Longitudinal and Lateral-Directional EOM in State Space Form 

The full 8x8 state space equation is shown below: X’=AX+BU. The longitudinal and 
lateral-directional EOM are decoupled in this analysis. The longitudinal EOM need not be solved 
to analyze this problem. The longitudinal EOM are included for completeness. 
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4. Determination of the Non-Dimensional Stability Derivatives 

The non-dimensional force and moment coefficients were determined in wind tunnel 
testing and updated based on actual flight data. The equations and coefficients are found in the 
Operational Aerodynamic Data Book Volumes I and III [Ref. 12]. 

a. Longitudinal Stability Derivatives 

The Taylor series first order expansion of aerodynamic coefficients is given by: 

Cd\ = Cdo + ('l)aCH + Cd*(&)\ 

Cl, = Clo + CLatn + Cl&(&) i 

Cm, = Cmo + ( mail 1 + C'//kv( A’) 1 

The elevons are at trail (0°) at touchdown [Ref 2]. Postlanding, the elevon is 
deflecting as required to balance the longitudinal moment equation and maintain a constant pitch 
attitude. However, elevon deflection is minimal until the start of derotation. At this time, the 
elevons ramp quickly downward as theta decreases at a rate of 1-2 deg/sec. Therefore, (5 e )i is 
negligible in the two point stance and set to 0 in these three equations. 

The coefficients [Ref. 1 2] are defined as: 

C(), — C( )BASIC ^ AC()high all ~^AC()real gas interaction ""^ACqbF ""("ACoeievator aileron AC( )plunging 

+AC ( )p lt ch rate "^ACqlG "^ACqGE AC( ) s ilts pod “t~AC( Vehicle 

where the subscript ( ) denotes lift force (L), drag force (D), or pitching moment (m) and: 
Q)basic = Basic, full-scale, freestream, rigid body force or moment coefficient. 

AC ()hlgh a i t = Change in force or moment coefficient due to high alt. effects (h>285,000 ft). 
AC ()rea i gas = Change in force or moment coefficient due to real gas effects. 

AC ()vl scous interaction = Change in force or moment coefficient due to viscous interaction (V* effects). 
AC () bf = Change in force or moment coefficient due to body flap deflection. 

ACoeievator aileron = Change in force or moment coefficient due to elevator and aileron deflection. 
AC( )p iun g ing = Change in force or moment coefficient due to rate of change of alpha. 

AC ( )p llch nuc = Change in force or moment coefficient due to pitch rate. 

AC () lg = Change in force or moment coefficient due to fully extended landing gear. 

AC()ge = Change in force or moment coefficient due to proximity of the ground. 

AC() S ii tspo d := Change in force or moment coefficient due to addition of silt pods. 
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AC()vehicie = Change in force or moment coefficient due to geometry differences between vehicles. 

AC 0 h, g hait and AC 0rca \ gas are not applicable below an altitude of 265,000 ft. and M 
< 6.0, respectively. A C ()visC ous interaction is not a factor since V® is zero. AC ()c icvaior/aiicnm is assumed to 
be zero since elevator/aileron deflection is considered negligible. AC()pi ung ing and AC() pitc h ra te are 
not a factor with the Orbiter on the ground in the two point stance since the pitch rate and angle of 
attack are essential constant until the initiation of derotation. This analysis will assume Orbiter 
Vehicle (OV) 103 which does not have a silts pod; therefore' AC( )s ii t s pod is not applicable. The 
above equation can be reduced to: 

C 0 i = C()BASIC + AC()BF“ h AC()LG + AC()GE ++ AC()vehicle 

Each of the terms in the above equation was obtained by linear interpolation of 
data extracted from a table at a = 10°, M=0.4 and a = 10°, M=0.25 to obtain data applicable for 
M=0.325 (353 fps) for both the lift and drag coefifiecients. Cia was determined by breaking each 
of C()basic> AC () bf, AC () lg, AC () ge, and AC {)V ehicie terms into Clo and Cl<x components and summing 
the contributions. C Da was determined in the same manner. 

Since U is considered constant in this analysis, C D and and C L are zero. 

C D and C L are zero since the pitch rate is zero. The lift and drag coefficients pertinent to this 
analysis are summarize in Table 3. 



Table 3. Lift and Drag Coefficients 



Lift and Drag Coefficients 


Value (per °) 


Cl i 


0.54055 


Cl<x 


0.05948 


Cd i 


0.16712 


Cdcc 


0.00992 



b. Lateral-Directional Stability Derivatives 

The Taylor series first order expansion of the lateral-directional aerodynamic 
coefficients are provided for completeness. The Taylor series first order expansion of 
aerodynamic coefficients is given by: 

Cyi = Cyo + Cyp/3\ + CySa(Sa) I + Cy&i^Sr) 1 



35 



Ch = Clo + Cl/l/3\ + Cl6o(Sa)\ + Cld(Sf)\ 

Cm = (no + ( nfij3l + Cn&(&) 1 + Cn&(&)i 

The lateral-directional bias coefficients are assumed to be zero for a symmetrical 
airframe: Cyo = Cio = Cno= 0. At touchdown, a software flag denoted ‘flat turn’ is set [Ref. 13] 
When this flag is set, lateral motion is accomplished without elevon motion, using rudder only. 
Therefore, aileron deflection is negligible in the two point stance and (&)i are set to zero in the 
above three equations. The coefficients [Ref. 1 2] are defined as: 

C()sidesUp +AC ()aileron ~^AC( (rudder ""^AC() ro u rate AC()y aw rate ^"^C()silts pod ~^AC() V ehicle 

where the subscript ( ) denotes side force (y), yawing moment (n), or rolling moment (1), and: 
Q)sidesiip = Basic, full-scale, ffeestream force or moment coefficient due to sideslip angle. 

AQ)aiieron = Change in force or moment coefficient due to aileron deflection. 

AC() m dder = Change in force or moment coefficient due to rudder deflection. 

AC ()ro u rate = Change in force or moment coefficient due to roll rate. 

AC()y awrat(: = Change in force or moment coefficient due to yaw rate. 

AC( )sl i, s ^ = Change in force or moment coefficient due to addition of silt pods. 

AC( ^vehicle = Change in force or moment coefficient due to geometry differences between vehicles. 
ACoaiieron may be dropped from this equation since aileron deflection is negligible. 

Each of the lateral-directional stability and control derivatives was obtained by 
linear interpolation of data extracted from a table at a = 10°, M=0.4 and a = 10°, M=0.25 to 
obtain data applicable for M=0.325 (353 fps). The side force, rolling moment, and yawing 
moment coefficients pertinent to this analysis are summarize in Table 4. 
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Table 4. Side Force, Rolling Moment, and Yawing Moment Coefficients 



Side Force, Rolling Moment, 
and Yawing Moment 
Coefficients 


Value (per °) 




-0.01705 




0 


Cy. 


0 




0.00138 


c >, 


-0.00205 


C 'p 


-0.3189 


C K 


0.1890 




0.0032 


c », 


0.0010 




0.0280 


C », 


-0.2925 




-0.00054 
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V. PROPOSED FCS MODIFICATIONS TO IMPROVE 



DIRECTIONAL CONTROL IN THE TWO POINT STANCE 

The rudder FCS susceptibility to PIO results from the slow surface rate, slow response 
time from rudder pedal input to rudder deflection, the power of the rudder, and the placement of 
the rudder with respect to the gear. Increasing the surface rate requires hardware modification. 
Analysis of hardware modifications are beyond the scope of this thesis. Transport delay might be 
decreased by shortening time delays in the FCS, linearizing the rudder pedal command shaper, 
and/or disabling Ny feedback. The power of the rudder might be reduced by decreasing the 
rudder pedal gain. 

Options to improve directional control may be categorized into expensive and 
inexpensive changes. Inexpensive options include those changes which require only software I- 
load or K-load change or a new piloting technique. Software I-loads and K-loads are constants 
which are not coded in the software and can be easily changed with no actual software code 
change. Expensive options include hardware changes, software changes which require a code 
change, and changes which require FCS recertification. 

Table 5 summarizes the problems with the current FCS, potential changes to the FCS to 
reduce susceptibility to PIO, the relative cost of the change, and whether the change is evaluated in 
this thesis. 



Table 5. Potential Changes to the FCS 



Problem 


Potential Fixes 


Relative Cost 


Evaluated 


Slow rudder rate 


Hardware modification 


Very expensive 


No 


Slow rudder response 


Decrease time delays 


Inexpensive 


Yes 




Linearize rudder command shaper 


Inexpensive 


Yes 




Disable Ny feedback 


Expensive 


Yes 


Powerful rudder 


Decrease pedal gain 


Inexpensive 


Yes 


Position of rudder 


Redesign Orbiter 


Very expensive 


No 
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A. 



LINEARIZE THE RUDDER PEDAL COMMAND SHAPER 



The rudder command shaper was described in Section III. A. 1 .b. Recall that the shaping 
schedule is parabolic and conforms to the following function: 

DRMS = (0.042 )DRMAN 2 + (0\2\)DRMAN . 

The shaper can be linearized by changing 0.042 to 0 and 0. 13 1. to 1.0. The above 
equation becomes: DRMS = DRMAN . This change provides a one-to-one response up to the 
limits of ± 22.5° which may result in a more predictable rudder response. 

B. REDUCE THE RUDDER PEDAL GAIN 

The rudder pedal gain was described in Section m.A.l.c. Reducing the gain will reduce 
the g-command authority on the rudder and may reduce over-controlling of the Orbiter. The 
current gain value is 0.08574 for a touchdown speed of 207 KEAS. The proposal is to limit the 
gain, GNYDRM to 0.04 at 683 fps TAS and below. This change will reduce the gain by 47% 
which may reduce the over-controlling tendency. 

C. DECREASE TIME LAGS 

The time lag between rudder pedal input and rudder deflection can be reduced by 
changing the Ny feedback time constant, described in Section III. A. 1 .f, from x = 0.8 s to x = 0.2 s. 
This change should result in quicker rudder response from the FCS which should reduce the 
susceptibility to PIO. 

D. DISABLE LATERAL ACCELERATION FEEDBACK 

Another option is to disable Ny feedback in combination with changing the Ny feedback 
time constant from t=0.8 s to x=0.2 s. This change would result in a more direct rudder FCS by 
allowing the rudder command to pass directly through the system with minimal delay. This 
change would require software code change and is relatively expensive. 
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VI. YAW FCS SIMULATION 



Directional control was evaluated by incorporating Orbiter lateral-directional control laws 
and longitudinal and lateral-directional EOM in MATLAB 5.0 and Simulink. The script program, 
'eom.m', defines the initial conditions, mass properties, stability and control derivarives, and 
longitudinal and lateral-directional EOM in state space form. The code for 'eom.m' is contained in 
Appendix A. The code for the Simulink simulation is contained in Appendix B. The plotting 
program, 'FCSplot.m', is contained in Appendix C. 

A. SIMULATION INITIAL CONDITIONS 

The case investigated is comparable to the STS-62 touchdown flight conditions. The 
commander input approximately +15° (right) mdder pedal deflection followed by -15° (left) 
rudder pedal deflection in the two point stance in an attempt to capture mnway centerline. 
Landing speed was approximately 207 KEAS with minimal crosswind (< 2 knots). Table 6 
contains the initial conditions used in the simulation. 



Table 6. Initial Conditions 



Parameter 


Initial Value 


Vehicle Weight 


228237 lbs 


Xcg 


1084.15 inches 


Ycg 


0 


Zcg 


370.65 inches 


Ixx 


8.58677e+05 slugs*!! 2 


lyy 


6.8699075e+06 slugs*!! 2 


Izz 


7.1334485e+06 slugs*!! 2 


Wing Surface Area 


2690 ft 2 


Chord Length 


39.56 ft 


U1 


207 KEAS (353 fps) 


Runway rolling \i 


0.03 


Sea Level Density 


0.002377 slugs/!! 3 


©, 


9° 


a 


10° 
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B. FCS CONFIGURATIONS EVALUATED 

Simulations were run to compare six FCS configurations as defined in Table 7. 



Table 7. FCS Configurations Evaluated 



Case Number 


FCS Configuration (Simulation Name) 


0 


Current configuration (baseline) 


1 


Disable Ny feedback (nony) 


2 


Decrease time lag t to 0.2 (tau) 


3 


Linearize rudder shaper (linmd) 


4 


Reduce rudder pedal gain (rudgain) 


5 


Cases 3,4,5 combined (done) 



C. SIMULATION 

The simulations were run in the time domain for 30 seconds using Dormand-Prince 
integrator. The simulated rudder pedal deflection included a rudder doublet of +15° rudder at 2 
seconds followed by -15° as depicted in Figure 19. The master Simulink block diagram is shown 
in Figures 20. 
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Figure 19. Rudder Pedal Deflection 
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Figure 20. Simulink Block Diagram of Rudder Doublet Effect on Directional Control 

The lateral-directional FCS is a subsystem of the entire system shown in Figure 20. The 
FCS is labeled TAW FCS’ in Figure 21. This subsystem contains the simulation code and block 
diagrams for the FCS. This block diagram and its internal block diagrams are shown in Figure 2 1 . 



43 






Figure 21. YAW FCS Block Diagram 
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VII. RESULTS 



Traces of 8^, 8 r , Ny, P, r, p, \\t and 4> versus time are shown in Figures 22-28. All six 
configurations are shown on each plot for comparison of the proposed changes. All Simulink 
traces shown in Figures 22-28 are in response to the 15° doublet. The Simulink model was 
validated by comparing the results with actual Shuttle flight data and Shuttle Entry Simulator 
(SES) simulation results. 

A. MODEL VALIDATION 

To establish validity of the Simulink model, actual STS-62 flight data and SES simulation 
results were compared to the Simulink results. The STS-62 flight data and SES results are 
contained in Appendix D. The STS-62 data is compared to the Simulink ’Baseline’ traces, since 
the onboard FCS is in the same configuration. The STS-62 CDR input a comparable 15° rudder 
pedal doublet in the two point stance. An SES evaluation, using the same initial conditions as this 
simulation and the STS-62 flight data, was run for comparison. Two SES runs were made. The 
SES ’Old FCS’ traces contain the same FCS configuration as the Simulink ’Baseline’ run. The SES 
'CF FCS' (Combined Fix) traces contain the same FCS changes as the Simulink 'Optimum fix'. 

Rudder deflection traces for all three data sets (Simulink, STS-62 flight data, and SES) 
show similar magnitude of nearly 15° deflection, similar response time, and deflection rate. Ny 
traces exhibit non-minimal phase response as indicated by the initial low magnitude acceleration 
in the negative direction. Lag times are similar and the magnitude of the Ny traces are 
comparable. A comparison of the sideslip angle traces indicate the same initial order-of- 
magnitude P resulting from the doublet. However, the Simulink trace shows an apparent long 
term directional instability which may be the reason why there is no P response in the opposite 
direction. The long term directional instability is reported to exist in the actual Orbiter in the two 
point stance. This instability is most likely not as strong as indicated in the Simulink model. 

Simulink yaw rate and roll rate traces also indicate the same magnitude and lag time as 
SES results and actual flight data. The SES ’Old FCS’ roll rate data is 'noisy' due to the gear 
model. 

This comparison of flight data and SES results with these Simulink traces results in a high 
level of confidence that this model is sufficient to analyze and compare the proposed FCS 
configurations for short term Orbiter directional response. 
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B. RUDDER DEFLECTION 



The rudder deflection plot in Figure 22 indicates that decreasing the time constant results 
in the quickest rudder response. Simply reducing x to 0, as in the 'No Ny feedback' option, yields 
a rudder response very similar to the pedal command. This option appears the best with respect to 
rudder surface control, yet the absence of Ny feedback generates other long term anomalies in 
yaw rate and siseslip stability. 

Reducing x to 0.2 with Ny feedback, as in the 'Shorter time constant' trace, yields an 
acceptable rudder surface response with the benefits of Ny feedback. This trace shows 
approximately one half second of time delay by the end of the four second doublet. This time 
delay could also be considered as -45° of phase lag. While 45° of phase lag is not desired, this 
option appears to be much better than some of the other options. The baseline configuration 
shows nearly two seconds of time delay before the rudder returns through the zero deflection 
position. This delay equates to nearly 180° of phase lag by the end of the maneuver. This 
response is obviously an unacceptable control response to a low frequency command. 

Linearizing the rudder pedal shaper, with no additional changes, would result in an 
increase in rudder deflection of 25% over the baseline. This larger rudder deflection for the 
relatively small rudder pedal command would aggravate the already difficult fine directional 
control. Linearizing the rudder with no other changes is not recommended. 

Reducing the rudder pedal gain, with no additional changes, is effective in decreasing the 
magnitude of the rudder surface response. However, the rise time in the deflection is too slow. 
The slow rise time distorts the predictability of the rudder response. This change does not effect 
the excessive phase lag. 

The optimum fix', which includes the x=0.2 time constant, linearized rudder shaper, and 
reduced rudder gain, result in the best shaped rudder response while retaining the Ny feedback. 
The deflection is shaped more like the rudder pedal input than all other Ny feedback cases. 
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Figure 22. Rudder Surface Deflection 



C. LATERAL ACCELERATION AT THE COMMANDER S STATION 

The lateral acceleration plot shown in Figure 23 shows that eliminating Ny feedback 
results in greater Ny gain and more time required to return to a zero Ny condition. This option 
displays the worst long term directional stability characteristics. 

Both the 'Shorter time constant’ and 'Linearized rudder' fixes also result in an increase in g 
by 25% over the baseline FCS. This increase in lateral acceleration control power could lead to 
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ty(g) 



directional oversensitivity. The 'Shorter time constant" fix results in the least phase lag. while 
'Linearized rudder’ yields the longest time delay of the Ny feedback options. 'Reduced rudder 
gain' results in an appropriate decrease in Ny gain. However, this change alone results in 
excessive time delay. 

The 'Optimum fix' has the least time lag and acceptable Ny gain characteristics. 

Lateral Areleration at Conmnender Station 
0.25 



0.2 

0.15 



01 



005 

0 



-0.05 
- 0.1 
-0.15 

0 5 10 15 2D 25 30 

Time (sec) 

Figure 23. Lateral Acceleration at the Commander's Station 
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D. SIDESLIP 



The sideslip plot shown in Figure 24 shows the total rudder control power in the classic 
sense of degrees of sideslip per degree of rudder pedal. Most modem aircraft exhibit rudder 
control power on the order of 1 :4, or one degree of sideslip for four degrees of rudder pedal. The 
'Baseline* Orbiter trace depicts rudder control power of 1:2.3. This apparent oversensitivity can 
lead to difficulty in precisely controlling sideslip. 

The 'No Ny feedback' option exhibits excessive rudder control power with an additional 
long term instability. This long term response is unacceptable. The change that optimizes the 
classic 1 :4 rudder control power ratio is the 'Optimum fix'. 

The sideslip plot also reveals the inherent long term directional instability of the Orbiter in 
the two point stance. This instability is due to the change in the center of rotation from the CG to 
the main wheels as discussed earlier. 
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Figure 24. Sideslip Angle 



E. YAW RATE AND HEADING 

Yaw rate and heading traces are shown in Figures 25 and 26. Heading was derived by 
integrating yaw rate. Long term directional instability is evident in both yaw rate and heading 
traces. The ’Shorter time constant' and 'Optimum fix' options indicate the least phase lag of all 
options with Ny feedback. 
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Figure 25. Yaw Rate 
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F. ROLL RATE AND ANGLE OF BANK 

Adverse roll is evident as shown in Figures 27 and 28. This response is undesirable 
because adverse roll contributes to the apparent reduced directional control. The 'Reduced rudder 
gain' results in the least adverse roll. However, this option results in significant phase lag. The 
"Optimum fix’ offers reduce adverse roll and the least phase lag. The long term positive dihedral 
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is evident in the long term positive roll rate. This roll rate is in response to the long term negative 
sideslip. 
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Figure 28. Angle of Bank 



G. NATURAL FREQUENCY AND DAMPING RATIO 

The so called 'directional natural frequency', gjd, and damping ratio, Co, were estimated 
from the plots after completion of the doublet. The damping ratio was estimated using the 
following equation: 
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= sin tan 



-i 



. D 

In — 

D, 



n 



Estimations from the plots indicate that tn D = 0.6 rad/sec and Q) = 0.45. As expected, the 
natural frequency and damping ratio are low. This low natural frequency and low damping ratio 
can result in the tendency of the FCS to be susceptible to PIO. 
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VIII. CONCLUSIONS 



Orbiter landings indicate both long term directional instability and the potential for pilot 
induced oscillations during landing and rollout before nosewheel touchdown with the baseline 
(current) FCS. The Orbiter’s FCS requires improvements to increase directional control in the 
two point stance. This thesis described the control deficiency, potential improvements to the FCS, 
and evaluated six FCS configurations to determine their effect on directional control. The 
evaluation was performed by modeling the Orbiter’ s postlanding lateral-directional control laws 
using MATLAB 5.0 and Simulink. The longitudinal equations of motion were include so that this 
model can more readily be adapted to evaluate longitudinal control problems. 

A. RECOMMENDATIONS 

Of the six FCS configurations evaluated, the option which provides the most 
improvement over the baseline configuration is the ’Optimum fix’. This fix includes reducing the 
time lag x from 0.8 sec to 0.2 sec, linearizing the rudder pedal shaper, and reducing the rudder 
gain. This change minimizes the apparent phase lag, minimizes the adverse roll, reduces 
directional oversensitivity, and provides the most traditional rudder control power. This option is 
relatively inexpensive since only software K-load and I-load changes are required. Software code 
and hardware modifications are not required. The data indicates that PIO tendencies would be 
reduced over the baseline system if this option were implemented 

This simulation was limited to analyzing the response of the FCS to a rudder doublet 
since the initial Orbiter directional difficulties were in response to a CDR doublet. The doublet 
input, in essence, is a frequency domain input. Therefore, frequency domain analysis would be 
appropriate. However, the validity of the Simulink model can more accurately be determined by 
comparing the time domain response traces of the individual states to actual flight data. Although 
this method was used as a diagnostic tool to study a frequency-concerned phenomenon, phase lag 
can be estimated from the traces. In addition, these changes can be comapared and evaluated with 
a pilot in the loop in time-based simulators such as the SES and Ames VMS. 

This author recommends further investigation of the 'Optimum fix' in the SES and Ames 
VMS. Frequency domain analysis would also be appropriate. In addition, the 
MATLAB/Simulink program used for this evaluation should be upgraded to include a gear model, 
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drag chute model, and wind model. The roll and pitch channels should be incorporated for 
completeness. 

B. USE OF COTS SOFTWARE FOR ENGINEERING EVALUATION 

This author concludes that commercially available engineering software packages such as 
MATLAB may be used to analysis proposed FCS changes prior to full up simulation in more 
expensive simulators. Future programs, such as X38, as well as Orbiter upgrade projects can 
benefit from the use of PC-based simulation. 
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APPENDIX A. MAI LAB CODE FOR SCRIPT PROGRAM 



This appendix contains the Matlab code for the script program, "eom.m". EOM.M 
defines the initial conditions, mass properties, stability and control derivatives, and longitudinal 



and lateral/directional EOM in state space form. 



^Initial conditions and EOM 
clear 
S = 2690 
C =39.56 
B=78 .057 
W=22 8237 
g=32 .174 
M=W / g ; 



^surface area of the wing in ft /v 2 

%chord length in ft 

%wing span in feet 

%weight in lbs 

%gravitational constant 

%mass in slugs 



Ixx=858677 ; 
Iyy=6869907.5; 
Izz=7133448 . 5; 



%Shift the moments of inertia from the eg to center of rotation 
% (between the main gear) 

I zmg=Izz+M* ( (1172.6-108 4.15) /12) A 2 

%b/h=0.05 at alpha = 10 

Mu=0.03; % runway rolling Mu 

Ul=353; %initial velocity in the x direction ft/sec 

rho= 0.002377; %sea level pressure 
theta=9*pi/180 %Pitch angle in radians 
alpha=10*pi/180 %Angle of attack in radians 
q=0 . 5*rho*Ul A 2 ^Dynamic pressure at sea level 

%Non-Dimensional Stability and Control Derivatives 
CLa=0. 05948 ; 

CDa=0. 009924 ; 

CL1=0. 54055; 

CD1=0. 16712; 

CLu=0; 

CDu=0; 



Cyb=-0 .01705; 

Cyp=0; 

Cyr=0; 

Cydeltar=0. 00138; 

Clb=-0. 00205; 

Clp=-0 .3189; 

Clr=0 . 1890 ; 

Cldeltar=0. 00032 ; 

Cnb=0 . 0010 ; 

Cnp=0 . 028 ; 

Cnr=-0 . 2925 ; 

Cndeltar=-0 . 00054 ; 

%Longitudinal Dimensional Stability and Control Derivatives 
Xu=-q*S* (CDu+2*CDl ) / (M*U1) 

Xa=-q*S* (CDa+2*CLl ) /M 
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Xq=0; 

Xdeltae=0; 

Zu=-q*S* (CLu+2*CLl ) / (M*U1) 

Za=-q*S* (CLa+2*CDl ) /M 
Z q=0; 

Zdeltae=0; 

Mu= 0 ; 

Ma=0; 

Mq=0 ; 

Mdeltae=0; 

%Lateral-directional Dimensional Stability and Control Derivatives 
Yb=q*S*Cyb/M 

Yp=q*S*B*Cyp/ (2*M*U1) %is zero 

Yr=q*S*B*Cyr/ (2*M*U1) %is zero 

Ydeltaa=0; 

Ydeltar=q*S*Cydeltar/M 

Lb=q*S*B*Clb/ Ixx 
Lp=q*S*B A 2*Clp/ (2* Ixx*Ul ) 

Lr=q*S *B A 2*Clr/ (2*Ixx*Ul) 

Ldeltaa=0; 

Ldeltar=q*S*B*Cldel tar/Ixx 

Nb= q * S * B * Cnb/ 1 zmg 
Np=q*S*B A 2*Cnp/ (2*Izmg*Ul) 

Nr=q*S*B A 2*Cnr/ (2*Izmg*Ul) 

Ndeltaa=0; 

Ndeltar=q*S*B*Cndeltar/Izmg 
%Definition of dx/dt= Ax+Bu 

%xdot=[udot alphadot qdot thetadot betadot pdot rdot phidot] 

%x=[u alpha q theta beta p r phi] 

%U= [ DELTAE DELTAA DELTAR] 

A = [Xu Xa Xq -g^cos (theta) 0000; 

Zu/Ul Za/Ul Zq/Ul+1 -g* sin ( theta) /U1 0000; 

Mu Ma Mq 0 0 0 0 0; 

0 0 1 0 0 0 0 0 ; 

0000 Yb/Ul Yp/Ul Yr/Ul-1 g+ cos ( theta) /U1 ; 

0 0 0 0 Lb Lp Lr 0; 

0 0 0 0 Nb Np Nr 0; 

000001 tan (theta) 0] 

B=[Xdeltae 0 0; 

Zdeltae/Ul 0 0; 

Mdeltae 0 0; 

0 0 0 ; 

0 Ydeltaa/Ul Ydeltar/Ul; 

0 Ldeltaa Ldeltar; 

0 Ndeltaa Ndeltar; 

0 0 0 ] 

% dy/dt = Cx+Du; C= identity matrix D=0 so that the output 
%of the state space is X 
C= [eye (8) ] 
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D= [ zeros (8,3) ] 

B= [Xdeltae 0 0; 
Zdeltae/Ul 0 0; 

Mdeltae 0 0; 

0 0 0 ; 

0 Ydel taa/Ul Ydeltar/Ul 
0 Ldeltaa Ldeltar; 

0 Ndeltaa Ndeltar; 

0 0 0 ] 
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APPENDIX B. MATLAB CODE FOR SIMULATION 



The code for the actual simulation, "baseline.m" is contained in this Appendix. The code 
is for the baseline system. 

Model { 



Name 


’’baseline" 


Version 


2.00 


SimParamPage 


Solver 


SampleTimeColors 


off 


WideVectorLines 


off 


PaperOrientation 


landscape 


StartTime 


"0.0" 


StopTime 


"30" 


Solver ode45 




RelTol M 01" 




AbsTol "le-6" 




Refine "1" 




MaxStep 


".05" 


InitialStep 


"auto" 


FixedStep 


"auto" 


MaxOrder 


5 


OutputOption 


RefineOutputTimes 


OutputTimes 


M [f 


LoadExtemallnput 


off 


Extemallnput 


”[t, u]" 


SaveTime 


off 


TimeSaveName 


"tout" 


SaveState 


off 


StateSaveName 


"xout" 


SaveOutput 


off 


OutputSaveName 


"yout" 


LoadlnitialState oflf 




InitialState 


"xlnitial" 


SaveFinalState off 




FinalStateName "xFinal" 




Limit MaxRows 


off 


MaxRows 


"1000" 


Decimation 


M j M 


AlgebraicLoopMsg 


warning 


MinStepSizeMsg 


warning 


UnconnectedlnputMsg 


warning 


UnconnectedOutputMsg 


warning 


UnconnectedLineMsg 


warning 


Consistencychecking 


off 


ZeroCross 


on 


BlockDefaults { 




Orientation 


right 


ForegroundColor 


black 
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BackgroundColor 


white 


Drop Shadow 


off 


NamePlacement 


normal 


FontName 


"Helvetica" 


Font Size 


10 


FontWeight 


normal 


FontAngle 


normal 


ShowName 


on 



} 

AnnotationDefaults { 



HorizontalAlignment 


center 


VerticalAlignment 




middle 


ForegroundColor 




black 


BackgroundColor 




white 


DropShadow 




off 


FontName 




"Helvetica" 


FontSize 




10 


FontWeight 




normal 


FontAngle 

} 

System { 




normal 


Name 


'baseline'* 


Location 




[0, 3, 644, 447] 


Open on 




ScreenColor 
Block { 




white 


BlockType 




Derivative 


Name 




"BETADOT" 


Position 

Orientation 

} 

Block { 


left 


[385, 320,415, 340] 


BlockType 




Clock 


Name 




"Clock" 


Position 




[30, 330, 50, 350] 


DeleteFcn 




"simclock BlockIsBeingDestroyed 


PostSaveFcn 


"simclock Save" 


Location 

\ 




[12,31, 119, 76] 


) 

Block { 






BlockType 




Demux 


Name 




"Demux" 


Position 




[370, 154,405,236] 


Outputs 




" g " 


Block { 






BlockType 




Mux 


Name 




"Mux" 


Position 




[310, 253, 335, 347] 


Orientation 

Inputs 


left 


"3 " 



64 



} 



Block { 






BlockType 




Mux 


Name 




"Muxl" 


Position 




[235, 178, 270, 212] 


Inputs 

\ 




”3" 


) 

Block { 






BlockType 




Scope 


Name 




"Ny" 


Position 

Orientation 


down 


[215, 395, 245,425] 


NamePlacement 




alternate 


Floating 




off 


Location 




[188, 235, 512, 474] 


Open 




off- 


Grid 


on 




TickLabels 

ZoomMode 


on 


on 


TimeRange 




"auto” 


YMin 




If ^ It 


YMax 




"5" 


SaveT oWorkspace 


off 


SaveName 




"ScopeData" 


LimitMaxRows 




on 


MaxRows 

Decimation 


h j M 


"5000" 


Samplelnput 


off 




SampleTime 

} 

Block { 


"0" 




BlockType 




Fen 


Name 




"Ny " 


Position 

Orientation 


left 


[240, 277, 275, 323] 


Expr 

} 

Block { 


"-(3 53 *u[3]+3 53 *u[2]+u[ 1 ]*60)/g' 






BlockType 




Integrator 


Name 




"PSI" 


Position 

ExtemalReset 


none 


[505, 380, 525, 400] 


InitialConditionSource 


internal 


InitialCondition 

LimitOutput 


off 


"0" 


UpperSaturationLimit 


"inf’ 


LowerSaturationLimit 


”-inf’ 


ShowSaturationPort 


off 


ShowStatePort 




off 


AbsoluteTolerance 


"auto" 



} 
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Block { 
BlockType 
Name 
Position 

ShowPortLabels 

MaskType 

MaskDescription 

MaskHelp 



MaskPrompt St ring 

MaskStyleString 

MaskVariables 

Masklnitialization 

MaskDisplay 



Subsystem 
" Pul seVnGenerator 1" 

[25, 35, 55, 65] 
on 

"Pulse Generator" 

"Pulse Generator" 

"Pulse generator which ensures pulse transitions " 
"are hit. Provides a vector of pulses when " 

"the height is entered as a vector. Unmask " 

"to see how it works." 

"Period (secs):|Duty cycle (% of period):|Amplitu" 
"de:|Start time:" 

"edit,edit,edit,edit" 

"period=@ 1 ;duty=@2;amplitude=@3;start=@4;" 
"\n" 

plot(0,0, 1 00, 1 00,[95,90,90,75,75,60, 60,45,45,30," 

"30, 15, 15, 5], [25, 25, 75,75, 25,25, 75,75, 25, 25, 75, 75 
",25,25]), \n 



•t 



Vn 



tf 



m 



MasklconFrame 

MasklconOpaque 

MasklconRotate 

MasklconUnits 

MaskValueString 

System { 

Name 

Location 

Open 

ScreenColor 
Block { 
BlockType 
Name 
Position 
Value 



on 

on 

off 

autoscale 

"1000|.2|15|5" 

"Pulse\nGenerator 1 " 
[205, 368, 544, 559] 
off 
white 

Constant 
"Constant 1" 

[40, 135, 60, 155] 
" 0 " 



Block { 
BlockType 
Name 
Position 
Amplitude 
Frequency 
Phase 

SampleTime 

} 

Block { 
BlockType 
Name 



Sin 

"Sine Wave" 

[35, 83, 65, 107] 

•• | H 

"2*pi./period" 

"-2*pi*(start./period + duty/200 - 1/4)" 



Step 

"Step" 
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Position 


[35, 27, 65, 53] 


Time 

Before 

After 


"start" 

"0" 

"amplitude" 


Block { 
BlockType 
Name 
Position 
Threshold 

} 

Block { 


Switch 

"Switch" 

[120, 65, 155, 125] 
"sin(pi*(0.5 - duty/100))" 


BlockType 

Name 

Position 


Outport 
"out 1 " 

[250, 88, 280, 102] 



Port "1" 

OutputWhenDisabled held 

InitialOutput "0" 



} 

Line { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 


"Sine Wave" 

1 

[5, 0; 20, 0] 

"Switch" 

2 


) 

Line { 
SrcBlock 
SrcPort 


"Switch" 

1 


DstBlock 

DstPort 


"outl" 

1 


} 

Line { 
SrcBlock 
SrcPort 


"Constant 1" 
1 


Points 

DstBlock 

DstPort 

} 

Line { 


[25, 0; 0, -30] 

"Switch" 

3 


SrcBlock 

SrcPort 

Points 

DstBlock 

DstPort 

} 

\ 


"Step" 

1 

[20, 0; 0, 35] 

"Switch" 

1 


) 

} 

Block } 
BlockType 
Name 


Subsystem 
" Pu lse\nGenerator2 " 
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Position 

ShowPortLabels 

MaskType 

MaskDescription 

MaskHelp 



MaskPromptString 

MaskStyleString 

MaskVariables 

Masklnitialization 

MaskDisplay 



[25, 95, 55, 125] 
on 

"Pulse Generator" 

"Pulse Generator" 

"Pulse generator which ensures pulse transitions " 
"are hit. Provides a vector of pulses when " 

"the height is entered as a vector. Unmask " 

"to see how it works " 

"Period (secs):|Duty cycle (% of period):|Amplitu" 
"de:|Start time:" 

"edit,edit, edit, edit" 

"period=@l;duty=@2;amplitude=@3;start=@4;" 

"\n" 

plot(0,0,100,100,[95,90,90,75, 75,60,60,45, 45,30," 

"30, 1 5, 1 5,5],[25,25,75,75,25,25,75, 75,25,25,75, 75 
M , 25,25]); \n 



fl 



II 

If 



\n 

fl 



MasklconFrame 

MasklconOpaque 

MasklconRotate 

MasklconUnits 

MaskValueString 

System { 

Name 

Location 

Open 

ScreenColor 
Block { 
BlockType 
Name 
Position 
Value 

} 

Block { 
BlockType 
Name 
Position 
Amplitude 
Frequency 
Phase 

SampleTime 



on 

on 

off 

autoscale 

"1000|.2|-15|7" 

"Pulse\nGenerator2" 
[205, 368, 544, 559] 
off 
white 

Constant 
"Constant 1" 

[40, 135, 60, 155] 
"0" 



Sin 

"Sine Wave" 

[35, 83, 65, 107] 
f'r 

"2*pi./period" 

"-2*pi*(start./period + duty/200 - 1/4)" 



Block { 
BlockType 
Name 
Position 
Time 
Before 



iiqm 



Step 

"Step" 

[35, 27, 65, 53] 



"start 



ii 
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After 



"amplitude" 



Block { 

BlockType 

Name 

Position 

Threshold 

} 

Block { 

BlockType 

Name 

Position 

Port 

Output WhenDisabled 
InitialOutput 

} 



Switch 

"Switch" 

[120, 65, 155, 125] 
"sin(pi*(0.5 - duty/ 100))' 



Outport 
"out 1 " 

[250, 88, 280, 102] 

H j H 

held 

" 0 " 



"Step" 

1 

0; 0,35] 
"Switch" 



Line { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 



[ 20 , 

1 



Line { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 

} 

Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 

Line { 
SrcBlock 
SrcPort 
Points 
DstBlock 
DstPort 



"Constant 1" 
1 

[25, 0; 0, -30] 

"Switch" 

3 



"Switch" 

1 

"out 1" 

1 



"Sine Wave" 



1 



[5, 0; 20, 0] 

"Switch" 



2 



} 

} 

Block { 
BlockType 
Name 
Position 
Orientation 

} 



Derivative 

"RDOT" 

[385, 260,415, 280] 
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Block { 

BlockType 

Name 

Position 

Floating 

Location 

Open 

Grid on 

TickLabels on 

ZoomMode 

TimeRange 

YMin 

YMax 

SaveToWorkspace 

SaveName 

LimitMaxRows 

MaxRows 

Decimation ”1" 

Sampielnput off 

SampleTime "0" 



Scope 

’’Scope” 

[395, 435,425,465] 
off 

[115, 431,439, 670] 
off 



on 

’’auto" 



off 

’’ScopeData’’ 

on 

”5000’’ 



Block { 

BlockType 

Name 

Position 

Floating 

Location 

Open 

Grid on 

TickLabels on 

ZoomMode 

TimeRange 

YMin 

YMax 

SaveToWorkspace 

SaveName 

LimitMaxRows 

MaxRows 

Decimation ’’1" 

Sampielnput off 

SampleTime ’’0" 

} 

Block { 

BlockType 

Name 

Position 

A 

B 

C 

D 



Scope 
’’ Scope 1’’ 

[390,380, 420,410] 
off 

[100, 161,424,400] 
off 



on 

’’auto’’ 



off 

’’ScopeData’’ 

on 



StateSpace 

"State-space" 



[295, 180,355,210] 
"A" 

"B" 

,’C” 

"D" 

X0 ”[0;0;0;0;0;0;0;0]” 



70 



} 

Block { 

BIockType 

Name 

Position 

Inputs 

} 

Block { 

BIockType 

Name 

Position 

VariableName 

Buffer 

Decimation " 1 " 

SampleTime "-I" 

} 

Block { 

BIockType 

Name 

Position 

VariableName 

Buffer 

Decimation " 1 " 

SampleTime H -1 M 

} 

Block { 

BIockType 

Name 

Position 

Orientation left 

VariableName 

Buffer 

Decimation "1" 

SampleTime "-1" 

} 

Block { 

BIockType 

Name 

Position 

VariableName 

Buffer 

Decimation "1" 

SampleTime "-1" 

} 

Block { 

BIockType 

Name 

Position 

VariableName 

Buffer 

Decimation " 1 " 



Sum 

"Sum" 

[100, 77, 130, 108] 



ToWorkspace 
"To Workspace" 
[95, 145, 155, 175] 
"doublet" 

"inf* 



ToWorkspace 
"To Workspace 1" 
[175, 250, 235,280] 
"dr" 

"inf 1 



ToWorkspace 
"To Workspace2" 
[95, 365, 155,395] 

"Ny" 

"inf’ 



ToWorkspace 
"To Workspace3" 
[540, 110, 600, 140] 
"beta" 

"inf’ 



ToWorkspace 
"To Workspace4" 
[570, 150, 630, 180] 

HpH 

"inf' 
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SampleTime "-l" 

} 

Block { 

BlockType 

Name 

Position 

VariableName 

Buffer 

Decimation "1" 

SampleTime M" 

} 

Block { 

BlockType 

Name 

Position 

VariableName 

Buffer 

Decimation M 1 " 

SampleTime M -l" 

} 

Block { 

BlockType 

Name 

Position 

VariableName 

Buffer 

Decimation " 1 M 

SampleTime ”-1” 

} 

Block { 

BlockType 

Name 

Position 

VariableName 

Buffer 

Decimation " 1 ” 

SampleTime H -l" 

} 

Block { 

BlockType 

Name 

Position 

ShowPortLabels 
System { 

Name 

Location 

Open 

ScreenColor 
Block { 
BlockType 
Name 



ToWorkspace 
"To Workspace5 M 
[585, 220, 645, 250] 

Mj.ll 

’’inf’ 



ToWorkspace 
"To Workspace6" 
[585, 335, 645, 365] 
’’psi” 

"inf' 



ToWorkspace 
’’To Workspace?" 
[70, 320, 130, 350] 
"time" 

"inf’ 



ToWorkspace 
"To Workspace8" 
[585, 280, 645, 310] 
"aob" 

"inf 1 



Subsystem 
"YAW FCS" 

[175, 179,205, 231] 
off 

"YAW FCS” 

[4, 44, 648, 440] 
off 
white 

Inport 
’’in 1” 
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Position 

Port 

Port Width 
SampleTime 

} 

Block { 
BlockType 
Name 
Position 
Port 

Port Width 
SampleTime 



[35, 180, 50, 200] 

H | It 

"-r 

m i K 



Inport 

"in_2" 

[100, 290, 120, 310] 

M 2 H 

"-r 



Block { 

BlockType 

Name 

Position 

Port 

Port Width 
SampleTime 

} 

Block { 

BlockType 
Name 
Position 
Orientation 
ShowPortLabels off 
System { 

Name 
Location 
Open 

ScreenColor 
Block { 

BlockType 
Name 
Position 
Orientation 
MATLABFcn 
OutputWidth 

} 



Inport 

"in_3" 

[290, 270,310, 290] 
”3" 

tl | M 
It | II 



Subsystem 

"COSALF" 

[320, 70,350, 120] 
down 



"COSALF" 

[318, 77, 607, 247] 



off 



down 



r 



white 

MATLABFcn 
M Cos( Alpha)” 

[57, 95, 133, 125] 

"cos" 



Block { 
BlockType 
Name 
Position 
Orientation 
Value 

} 



Constant 

"alpha" 

[58, 50, 132, 65] 

down 

’T0*pi/180" 



Block { 
BlockType 
Name 
Position 



Outport 

"out_P 

[85, 150, 105, 170] 
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Orientation down 

Port T 

Output WhenDisabled held 

InitialOutput "0" 



Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 



"Cos( Alpha)" 
1 

"out _1 M 
1 



Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 



"alpha" 

1 

"Cos( Alpha)" 
1 



} 

Block { 
BlockType 
Name 
Position 
Value 



Constant 

"Constant" 

[60, 125, 80, 145] 
"04" 



Block { 

BlockType 

Name 

Position 

ShowPortLabels off 
System { 

Name 
Location 
Open 

ScreenColor 
Block { 

BlockType 
Name 
Position 
Port 

Port Width 
SampleTime 



Subsystem 

"DRERR" 

[335, 255, 365, 305] 



"DRERR" 

[180, 294, 476, 476] 



off 



h | ii 
H j ii 



white 

lnport 

"inj" 

[15, 75, 35, 95] 

m j •• 



Block { 
BlockType 
Name 
Position 
Numerator 
Denominator 
SampleTime 

} 



DiscreteFilter 
" Y AWRATEBBF" 

[60, 56, 245, 114] 

"[.7718479 -.3633895 .4980655]" 
"[1 -.3633895 .2699135]" 
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Block { 

BlockType 

Name 

Position 

Port "P 

OutputWhenDisabled 
InitialOutput "0" 

} 

Line { 

SrcBlock 

SrcPort 

DstBlock 

DstPort 

} 

Line { 

SrcBlock 

SrcPort 

DstBlock 

DstPort 

} 

} 

} 

Block { 

BlockType 

Name 

Position 

Value 

> 

Block { 

BlockType 

Name 

Position 

Inputs "2" 

} 

Block { 

BlockType 

Name 

Position 

Inputs " 2 " 

} 

Block { 

BlockType 

Name 

Position 

Orientation 

Value 

} 

Block { 

BlockType 

Name 

Position 



Outport 

"outj* 

[265, 75, 285, 95] 

held 



" YAW_RATE_BBF" 
1 

"out_r 

i 



"in 1 " 

1 

M Y AW_RATE JBBF" 
1 



Constant 

M DRT M 

[135, 245, 150, 265] 
” 0 ” 



Product 

"GAIN" 

[135, 175, 160, 195] 



Product 

"GAIN" 

[355, 185, 380, 205] 



Constant 

"GDRC" 

[475, 300,515, 330] 
up 

"6.325" 



Constant 

"GERYFWDS" 

[485, 105, 505, 125] 
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Orientation down 

Value ".5" 



Block { 
BlockType 
Name 
Position 
Orientation 
ShowPortLabels 
System { 

Name 

Location 

Open 

ScreenColor 
Block { 
BlockType 
Name 
Position 
Orientation 
Input Values 
OutputValues 



Subsystem 
"GNYDRM_COMP" 
[105, 65, 135, 115] 
down 



off 



"GNYDRMCOMP" 

[355, 132, 620, 292] 
off 

white 

Lookup 

M GNY DRM Comp” 
[75, 43, 150, 107] 
left 

H [0 250 900 10000]" 

"[ 0.1 0.1 0.01 0.01 ] M 



Block { 

BlockType 

Name 

Position 

Orientation left 

Value 

} 

Block { 

BlockType 

Name 

Position 

Orientation left 

Port H r 

OutputWhenDisabled 
InitialOutput "0" 



Constant 

M TAS" 

[170, 75, 190, 95] 
"353" 



Outport 
"GNYDRM" 
[15, 70,35, 90] 



held 



Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 



"GNY DRM Comp” 
1 

"GNYDRM" 

1 



} 

Line { 
SrcBlock 
SrcPort 
DstBlock 
DstPort 

} 



"TAS" 

1 

"GNY DRM Comp” 
I 
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Block { 

BlockType Subsystem 

Name "NYCOMP" 

Position [295, 170, 325, 220] 

ShowPortLabels off 
System { 

Name "NYCOMP" 

Location [227, 109, 594, 284] 



Open off 

ScreenCoIor 
Block { 

BlockType 

Name 

Position 

Port T 

PortWidth 

SampleTime "-1" 

} 

Block { 

BlockType 

Name 

Position 

Gain "5" 

} 

Block { 

BlockType 

Name 

Position 

Numerator 

Denominator "[04 

} 

Block { 

BlockType 

Name 

Position 

Port "1" 
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APPENDIX C. MATLAB CODE FOR PLOT PROGRAM 



plot( time, doublet, 'b') 
title('Rudder Pedal Deflection') 
xlabel(’Time (sec)') 
ylabel('dpedal (deg)’) 
figure 

plot(time,dr, , - , ,time,drl, , ~ ' ',time,dr2,':',time,dr3,'x',time,dr4,\',time,dr5,'-.^ 
grid 

title('Rudder Deflection’) 
xlabel('Time (sec)’) 
ylabel('dr (deg)') 

legend('Baseline',’No Ny feedback','Shorter time constant', 'Linearized rudder', Heduced rudder 

gain','Optimum fix’, 1) 

figure 

plot(time,Ny,'-',time,Nyl,'— ',time,Ny2,':',time,Ny3,'x',time,Ny4,’.',time,Ny5,'-.') 
grid 

title('Lateral Acceleration at Commander Station') 
xlabel('Time (sec)') 
ylabelCNy (g)') 

legend( r Baseline', r No Ny feedback','Shorter time constant'^inearized rudder', 'Reduced rudder 

gain', 'Optimum fix’, 1 ) 

figure 

plot(time,beta* 1 8 0/pi,'-', time, beta 1 * 1 80/pi,'— 

',time,beta2* 1 80/pi,':', time, beta3 * 1 80/pi, 'x',time,beta4* 1 80/pi,'.', time,beta5* 1 80/pi,'-.’) 
grid 

title('Sideslip') 
xlabeI('Time (sec)') 
ylabelCBeta (deg)') 

legend( , Baseline','No Ny feedback','Shorter time constant', 'Linearized rudder','Reduced rudder 

gain', 'Optimum fix', 4) 

figure 

p!ot(time,r* 1 80/pi,'-', time,r 1 * 1 80/pi,'— 

', time, r2* 180/pi, ':',time,r3*l 80/pi, 'x’,time,r4* 180/pi, ’.',time,r5* 180/pi,'-.') 
grid 

title(’Yaw Rate') 
xlabel('Time (sec)') 
ylabel(’r (deg/sec)') 

legend('Baseline','No Ny feedback','Shorter time constant^ldnearized rudder’,’Reduced rudder 

gain',’Optimum fix',1) 

figure 

plot(time,p* 1 80/pi, '-',time,p 1 * 1 80/pi,’-- 

',time,p2* 1 80/pi, ':',time,p3 * 1 80/pi,'x’,time,p4* 1 80/pi,'.',time,p5* 1 80/pi,'-.') 
grid 

title('Roll Rate') 
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xlabel('Time (sec)') 
ylabel(’p (deg/sec)’) 

legend( , Baseline', , No Ny feedback','Shorter time constantVLinearized rudder','Reduced rudder 

gainVOptimum fix',1) 

figure 

piot(time,psi* 1 80/pi,'-',time,psi 1 * 1 80/pi,'~ 

',time, psi2M 80/pi,':', time, psi3*180/pi,’x',time,psi4*l 80/pi,'. ',time,psi5* 180/pi,'- ') 
grid 

title('Heading') 
xIabel('Time (sec)') 
ylabel('psi (deg)') 

legend('Baseline', r No Ny feedback','Shorter time constantVLinearized rudder','Reduced rudder 

gainVOptimum fix',1) 

figure 

plot(time,aob* 1 80/pi, '-',time,aobl * 1 80/pi,'— 

',time,aob2* 180/pi,':', time,aob3*180/pi,'x',time,aob4* 180/pi,'. ’,time,aob5* 180/pi,'-.') 
grid 

title(' Angle of Bank') 
xlabel('Time (sec)') 
ylabel('phi (deg)’) 

legend('Baseline','No Ny feedback’,'Shorter time constant', 'Linearized rudder'/Reduced rudder 
gain',’Optimum fix', 4) 
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APPENDIX D. STS-62 FLIGHT DATA AND SES SIMULATION 



RESULTS 

The first four plots contain the STS-62 data. The remaining six plots are the SES results. 
For all cases. Main Gear Touch Down (MGTD) occurs at 80.3 seconds with derotation at 91 
seconds. The X scale reflects time from the Approach and Land (A/L) transition, which normally 
occurs around 1 0,000 feet. 
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